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Abstract 

It is known that there are three maximally entangled states |<I>i) = (|0000) + |llll))/\/2, |<I> 2 ) = 
(A2|llll) + |1000) + |0100) + |0010) + |0001))/A6, and l^a) = (|1111) + |1100) + |0010) + |0001))/2 
in four-qubit system. It is also known that there are three independent measures (j = 1, 2,3) 
for true four-way quantum entanglement in the same system. In this paper we compute and 
their corresponding linear monotones for three rank-two mixed states pj = p|<hj)(<hj|-|-(l — 
p)|W 4 )(W 4 |, where IW 4 ) = (|0111)-|-|1011)-|-|1101)-|-|1110))/2. We discuss the possible applications 
of our results briefly. 
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I. INTRODUCTION 


Recently, much attention is being paid to quantum information theory (QIT) and quan¬ 
tum technology (QT) 1|. Most important notion in QIT and QT is a quantum correlation, 


which is usually termed by entanglement [2[ of given quantum states. As shown for last 
two decades it plays a central role in quantum teleportation superdense codingjd], quan¬ 
tum cloningj^, and quantum cryptography j^, 7|. It is also quantum entanglement, which 
makes the quantum computer^ outperform the classical one|9|. Thus, it is very important 


to understand how to quantify and how to characterize the entanglement. 


A. entanglement measures 


For bipartite quantum s yste m many entanglement measures were constructed before such 


as distillable entanglement llOl 


, entanglement of formation (EOF)[iy, and relative entropy 


of entanglement IREEI jlll. Il2|. 

The distillable entanglement is dehned to quantify how many maximally entangled states 
can be constructed from the copies of the given quantum state in the asymptotic region. 
Thus, in order to compute the distillable entanglement we should hnd the optimal purihca- 
tion (or distillation) protocol. If, for example, the optimal protocol generates n maximally 
entangled states from m copies of the quantum state p, the distillation entanglement for p 
is given by 

T) 

D{p) = lim -. (1.1) 

m^oo m 

Although the distillable entanglement is well-dehned, its analytical calculation is very dif- 
hcult because it is highly non-trivial task to hnd the optimal purihcation protocol except 
very rare cases [l^. 


REE of a given quantum state p is dehned as 


Er{p) = min^(p||cT), 

aev 


( 1 . 2 ) 


where "D is a set of separable states and S'(p||cr) is a quantum relative entropy; that is 
S'(p||(t) = tr(plnp — pin a). It is known that Er{p) is an upper bound of the distillable 


^ The current status of quantum computer technology was reviewed in Ref. [8| . 
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entanglement. However, for REE it is also highly non-trivial task to hnd the closest separable 
state a of the given quantum state p. Still, therefore, we do not know how to compute REE 
analytically even in the two-qubit system except rare cases jl^. 

The EOF for bipartite pure states is dehned as a von Neumann entropy of each party, 
which is derived by tracing out other party. For mixed state it is dehned via a convex-roof 
method Q,Q; 

= min^PjEi7’('0j), (1.3) 

j 

where minimum is taken over all possible pure state decompositions, i.e. p = 

with 0 < Pj < 1. The decomposition which minimizes is called the optimal 

decomposition. For two-qubit system, EOF is expressed as|l 6 | 


Ef{C) = h 


(1.4) 


where h{x) is a binary entropy function h{x) = —x\nx — (1 — x) ln(l — x) and C is called 
the concurrence. For two-qubit pure state \'ip) = 'ipij\ij) with {i,j = 0,1), (7 is given by 




(1.5) 


where the Einstein convention is understood and is an antisymmetric tensor. For two- 
qubit mixed state p the concurrence C{p) can be computed hy C = max(Ai — A 2 — A 3 — A 4 , 0), 
where {A^, A 3 , A 4 } are eigenvalues of p{ay ® o'y)p*{ay ® (Ty) with decreasing order. Thus, 

one can compute the EOF for all two-qubit states in principle. 

B. Classification of Entanglement 


Although quantihcation of the entanglement is important, it is equally important to 
classify the entanglement, i.e., to classify the quantum states into the different type of 
entanglement. The most popular classihcation scheme is a classihcation through a stochas¬ 
tic local operation and classical communication (SLOCC) 17|. If \'ijj) and |0) are in same 
SLOCC class, this means that l"^) and |0) can be used to implement same task of quantum 
information process although the probability of success for this task is different. Mathe¬ 
matically, if two n-party states \'ip) and |0) are in the same SLOCC class, they are related 
to each other by |' 0 ) = Ai ® y 42 ® ® A„| 0 ) with {Aj} being arbitrary invertible lo- 
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cal operators^. Moreover, it is more useful to restrict ourselves to SLOCC transformation 
where all {Aj} belong to SL(2, C), the group of 2 x 2 complex matrices having determinant 
equal to 1. In the three-qubit pure-state system it was shown [l^ that there are six differ¬ 
ent SLOCC classes, fully-separable, three bi-separable, W, and Greenberger-Horne-Zeilinger 
(GHZ) classes. Subsequently, the classification was extended to the three-qubit mixed-state 


system 


19|. 


The SLOCC transformation enables us to construct the entanglement measures for the 
multipartite states. As Ref. 20(] showed, any linearly homogeneous positive function of a 
pure state that is invariant under determinant 1 SLOCC operations is an entanglement 
monotone. One can show that the concurrence C in Eq. fll.Sp is such an entanglement 
monotone as follows. Let I'ljj) = 'ipij\ij) with i,j = 0, 1. Then, 1-0) = (A ® i?)|'0) = V'plo), 
where ifjij = xf^apAiaBj^. Using eijMiaMjp = (detM)ea /3 for arbitrary matrix M, it is easy to 
show = (cletA)(deti?)ej^j20ii2'l/^iiii'l/^*2i2) which implies that C is invariant 

under determinant 1 SLOCC operations. 

The theorem in Ref. 20|, i.e. a linearly homogeneous positive function that remains in¬ 
variant under determinant 1 SLOCC operation is an entanglement monotone, can be applied 
to the three-qubit system. If \ f)) = 'ipijkfjk), the invariant monotone is 

1/2 


o = 




( 1 . 6 ) 


This is exactly the same with a square root of the residual entanglement^ introduced i 


m 


Ref. 


2l| . The three-tangle fll.6p has following properties. If lip) is a fully-separable or a 


partially-separable state, its three-tangle completely vanishes. Thus, measures the true 
three-way entanglement. It also gives r 3 (GHZ 3 ) = 1 and TsiW^) = 0 to the three-way 
entangled states, where 


|GHZ3)^(|000) + |111)) 


|W3) = ^(1001) +1010)+ I100)), (1,7) 


For mixed state quantihcation of the entanglement is usually defined via a convex-roof 
method 3, Although the concurrence for an arbitrary two-qubit mixed state can be, in 


principle, computed following the procedure introduced in Ref. 16j, still we do not know how 


^ For complete proof on the connection between SLOCC and local operations see Appendix A of Ref. 181. 
^ In this paper we will call T 3 three-tangle and r| residual entanglement. 


4 

















to compute the three-tangle (or residual entanglement) for an arbitrary three-qubit mixed 
state. However, the residual entanglement for several special mixtures were computed in 
Ref. 22|]. More recently, the three-tangle for all GHZ-symmetric states 23(] was computed 


analytically 


Mi¬ 


lt is also possible to construct the SLOCC-invariant monotones in the higher-qubit sys¬ 
tems. In the higher-qubit systems, however, there are many independent monotones, because 
the number of independent SLOCC-invariant monotones is equal to the degrees of freedom 
of pure quantum state minus the degrees of freedom induced by the determinant 1 SLOCC 
operations. For example, there are 2(2” — 1) — 6 n independent monotones in n-qubit system. 
Thus, in four-qubit system there are six invariant monotones. Among them, it was shown 
in Ref. 2^ by making use of the antilinearity 1^ that there are following three independent 
monotones which measure the true four-way entanglement: 


• (o-^o-2(Tx(T2) • (aaO-'^cr^cra) 

^2^ = i.O-^aya2a2) • {<J^<J2<J\<J2) • (o-20-‘'(72(Tt) • {a2(J2(J^(J^) ( 1 . 8 ) 

•^3^^ = • (cr^o-^'aacra) • (o-pCTscr^aa) • (cr^’crscr^crs) • (a^craCTsaA) • (a^'cracraa^). 


where ai = I 2 , (Ti = cr^, a 2 = CTy, and the Einstein convention is introduced 

with a metric = diag{—1,1, 0,1}. Furthermore, it was shown in Ref. 2^ that there are 
following three maximally entangled states in four-qubit system: 


4>i> = ^(10000)+ |1111» 

■i>2) = -= f\/2|llll) + |1000> + |0100> + 10010) + |0001) 

Vb 

$ 3 ) = + |1100) + lOOlO) + lOOOl)). 



Tr(4) -77(4) -77(4) 

''l ''2 -^3 

l^l) 

1 1 I 

1 ^ 2 ) 

1 0 0 

1 ^ 3 ) 

0 0 1 

IW 4 ) 

0 0 0 


Table \ and of the maximally entangled and W 4 states. 


(1.9) 
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( 1 . 10 ) 


The measures J^ 2 ^\ and of |<hi), |<h 2 ), I'hs), and 

IW4) = ^(lOlll) + |1011) + |1101) + |1110)) 

are summarized in Table I. As Table I shows, |<hi) is detected by all measures while |<h 2 ) (or 
|<h 3 )) is detected by only (or As three-qubit system, IW 4 ) is not detected by all 

measures. 


C. Physical Motivations 


As states earlier, W and GHZ classes represent the true 3-way entanglement in three- 
qubit system. However, the three-tangle and residual entanglement r| cannot detect 
the entanglement of W class, but yield a maximal value to GHZ class. Then, it is natural 
to ask how much entanglement is 
p|GHZ 3)(GHZ3|+(1 -p)|W3)(W3 
whose residual entanglement is 


5 detected by rs and r| for the rank-2 mixture p{v) = 
|. This was explored in the first reference of Ref. j22| . 


where 


T'|(P(P)) 


0 

' 9i{p) 


for 0 < p < po 

for po < p < Pi 


gn{p) for Pi < p < 1 


( 1 . 11 ) 


9i{p) =P^ - ^ Vp(1 - Pf f///(p) = 1 - (1 -P) 0 + (1-12) 


Thus, one can say that the the influence of W class is dominant at 0 < p < po while influence 
of GHZ class is dominant at pi < p < 1. In the intermediate region Po < P < Pi two classes 
seem to compete with each other. For three-tangle similar method can be applied and the 
result is 


^3(Pp) 


0 for 0 < p < Po 
fzfr for Po < p < 1 . 


(1.13) 


The expression of r 3 is much simplier than that of r|. It is mainly due to the fact that rs is 
a linear invariant under the SLOCC transformation. 
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SLOCC classiheation 

representative states 

Gabcd 

*^2), h^s) 

Lab^ 

IW4) 

Labc2 

|0000) 

-^02 b2 

|0110)+ |0011) 

■^“20301 

|0011)+ |0101)+ |0110) 

■^O 30 jO 30 j 

|0000)+ |0111) 

-^ 7 ®I 

|0000)+ |1011)+ |1101)+ |1110) 


|0001)+ |0110)+ |1000) 

-^ 5®3 

|0000)+ |0101)+ |1000)+ |1110) 


Table II: SLOCC classification of four-qubit system. 


One can ask same question to four-qubit system by choosing a rank -2 state. However, 
situation is much more complicated because the four-qubit system has nine different SLOCC 
classes 27 (]. The nine classes and their representative states are summarized in Table II. Thus, 


there are too many combination to choose the rank -2 states. Motivated by the three-qubit 
case we choose two classes Gated and Lat^ and construct the rank-2 state pj = p|<hj)(d>j|-|-(l — 
p)|W4)(W4| (j = 1 , 2 , 3 ). The purpose of this paper is to compute and (j = 1 , 2 , 3 ), 
where is a linear entanglement monotone dehned as 












( 1 . 14 ) 


We will show that as in the three-qubit case the influence of Lat^ class is strong at 0 < p < poi 
where po is dependent on |<hj) and is larger than the corresponding 3 -qubit value 0.6269 for 
most cases. Of course, with increasing p the influence of Gated becomes stronger gradually. 

The paper is organized as follows. In sections II we derive the entanglement of pi, p2, 
and ps analytically. We also derive the optimal decompositions explicitly for each range 
in p. To check the correctness of our results we use the criterion discussed in Ref.| 28 |, i.e. 
entanglement should be a convex hull of the minimum of the characteristic curves. In section 
III we discuss the possible applications of our results. In the same section a brief conclusion 
is given. 
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II. ENTANGLEMENT OF pj (j = 1,2,3) 


In this section we will compnte the entanglement of pj = p|<I)j)(<I)j| + (l — p)|W 4 )(W 4 |. 
Before explicit calculation it is convenient to discuss the general method of our calculation 
briefly. First, we define a pure state 

\Zj{p,ip)) = -e*^y/r^|W4). (2.1) 

Since \Zj{p, ip)) is a pure state, one can compute {i = 1, 2, 3) of it by making use of Eq. 
(II. 8 p . As shown below, has a nontrivial zero at p = po for most cases. This is due to 
the fact that cannot detect the entanglement of IW 4 ). Exploiting this fact we construct 
the optimal decomposition, which yields = 0 at 0 < p < po. At po < p < 1 region 

we conjecture the optimal decomposition and corresponding J^l‘^\pj) by making use of the 
continuity of entanglement with respect to p and convex condition. Same procedure can 
be applied to the computation of Qf'\pj) {i,j = 1,2,3). Finally, we adopt a numerical 
method, which gurantees the correctness of our guess. 

A. Case pi 

In this subsection we will compute the entanglement of pi = p|<I>i)($i| + (l —p)|W 4 )(W 4 |. 
One can show 

[Zi(p, p)] = p|p2 - 3(1 - 

) [Zi(p, p)] = p >2 _ 4(1 _ ( 2 . 2 ) 

[Zi(p,p)] = ^, 

where \Zi{p, ip)) is defined in Eq. fl2.ip with j = 1. 

1. (pi) and (pi) 

From Eq. fl2.2p one can show that [Zi{Pi p)] has a nontrivial zero (<p = 0) 

x/3 

x/3 + 1 
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0.634. 


(2.3) 









The existence of finite po guarantees that {pi) should vanish at 0 < p < po- p = po 
this fact can be verified because we have the optimal decomposition 


Pi(Po) = 


1^1 {Po, 0))(2'i (po, 0) \ + \Zi ^po, —j)(Zi ^po, — 

+ |Zi(po,ir)){Zi(p(i,ir)| + |Zi {lHi,^))(Zi fpo, y 


(2.4) 


At the region 0 < p < po, {pi) should vanish too because one can find the following 
optimal decomposition 


P1(P) = —Pl(Po) + ( 1 - — ) |W4)(W4| 
Po \ Po, 


(2.5) 


Combining these facts, one can conclude that (pi) = 0 at 0 < p < po- 

Next, we consider the po < p < 1 region. Eq. fl2.4p at p = po strongly suggests that the 
optimal decomposition at this region is 


jZi {p,0)}(Zi (p,0) | + |Zi (^p,^'j}(Zi (^p,|) I 
+ jZi {p,7r)}(Zi{p,7r)j + jZi ^p, y^)(Zi fp,Y 


Pi(p) = ^ 


If Eq. Il2.till is a correct optimal decomposition in this region, J-i (pi) reduces to 

d"’ (Pi) = p(6p - 2p" - 3). 


(2.6) 


(2.7) 


Since the right-hand side of Eq. 
right. In conclusion, we can write 


is convex, our conjecture (Eq. fl2.6p j seems to be 


d*' (pi) = ^(p - Po)p(6p - 2p^ - 3), 


(2.8) 


where 0{x) is a step function defined as 


0{x) 


1 X > 0 

0 a; < 0. 


(2.9) 


However, if our choice Eq. 02.61) is incorrect, Eq. 02. 8 p is merely an upper bound of 
(pi). Thus, we need to prove that Eq. 02.81) is really optimal value. To prove this we 
should examine, in principle, all possible decompositions of pi, i.e. pi = and 

minimize the corresponding value However, it is impossible because pi has 

inhnite number of decomposition. 
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In order to escape this difficulty to some extent one may rely on Caratheodory’s theorem 


for convex hull 29(], which states that for four-qubit rank-2 states hve vector decomposition 
is sufficient to minimize J^^\pi). Thus, we need to investigate decompositions with 2, 3, 4, 


or 5 vectors. This method was used in the first reference of Ref. 


22j to minimize the residual 


entanglement of three-qubit rank-2 mixture. Still, however, it is difficult, at least for us, to 
parametrize all decompositions with first few vectors. 

In this paper, therefore, we will adopt the alternative numerical method presented in 


Ref. 


28[. We plot the p-dependence of [Zi(p, p)] for various p (See solid 


ines of Fig. 


28| showed. 


1(a)). These curves have been referred as the characteristic curves. As Ref. 

(Pi) is a convex hull of the minimum of the characteristic curves. Fig. 1(a) indicates 
that Eq. fl2.8p (thick dashed line) is really the convex characteristic curve, which implies that 
Eq. fl2.8l) is really optimal. This method was also used in the third reference of Ref. 2^ to 
minimize the residual entanglement of three-qubit rank-3 mixture. 

Now, let us consider (pi). It is easy to show that (pi) vanishes at 0 < p < 
Po due to the optimal decomposition Eq. (12.bh . If one chooses Eq. (12.6h as an optimal 
decomposition at po < p < 1, the resulting (pi) is not convex in the full range. Thus, 
we should adopt a technique introduced in Ref. 


IS 


Pi(p) = 


P-P0 | 

1 -Po 


<hi)(<hi|4 


22j. In this case the optimal decomposition 
1 


p 


Po 


Pi(Po), 


( 2 . 10 ) 


which results in (pi) = (p — Po)/(l — Po)- Combining all these facts, one can conclude 


d‘‘’(P.) = »(p-Po)^. (2.11) 

1 -Po 

To confirm that Eq. fl2.1ip is correct, we plot the characteristic curves [Zi(p, p)] for 
various p as solid lines and Eq. fl2.1ip as thick dashed line in Fig. 1(c). This hgure shows 
that Eq. fl2.1ip is convex hull of the minimum of the characteristic curves, which strongly 
supports the validity of Eq. fl2.11l) . 


2. (pi) and (pi) 

From Eq. fl2.2p one can notice that [-2^1 (p, p)] has a nontrivial zero (p = 0) 

= 0.667. (2.12) 

o 
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FIG. 1: (Color online) Plot of the p dependence of (a) [Zi{p,(p)]^ (b) F^^ [Zi{p,(p)], and (c) 

[Zi(p, (/?)] in Eq. (|2.2I) . We have chosen p from 0 to 27r as an interval 0.1. The thick dashed 
lines correspond to F^^ (pi) in Eq. (|2.8p . F^f^ {pi) in Eq. (I2.13D and {pi) in Eq. (12.IIP . These 
figures indicate that Eq. (12.8p . Eq. (|2.13P . and Eq. (12.lip are convex hull of the minimum of the 
characteristic curves. 

Thus, Eq. (12. 5 p and Eq. (12.6p with po = 2/3 can be the optimal decompositions for F^'^ (pi) 
at 0 < p < po and Po < P < 1, respectively. Then, the resulting F^'^ (pi) becomes 

(Pi) = - Po)p>^ - 4(1 - p)^]. (2.13) 

In order to confirm that our result (I2.13p is correct, we plot the characteristic curves for 
various p (solid lines) and Eq. (I2.13P (thick dashed line) in Fig. 1 (b). As Fig. 1(b) exhibits, 
our result (I2.13P is convex hull of the minimum of the characteristic curves, which strongly 
supports that Eq. (I2.13|) is really optimal one. 

Similarly, (pi) becomes Eq. (12.lip with changing only po to 2/3. The corresponding 
optimal decompositions are Eq. (12.51) at 0 < p < po and Eq. (I2.inp at po < p < 1, 
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respectively. Of course, we have to change po to 2/3. 


3. (pi) and (pi) 

Eq. fl2.2p shows that [Zi(p, p)] doe not have nontrivial zero. In addition, it is 
independent of the phase angle (p. This fact may indicate that there are inhnite number of 
optimal decompositions for (pi). The simplest one is 

Pi{p) = ^l^i(P,0))(Zi(p,0)| + ^|Zi(p,7r))(Zi(p,7r)|, (2.14) 

which gives {pi) = p®/2. If one chooses Eq. fl2.14p as an optimal decomposition for 
^ 3 ^^ (pi), it generates ^g"^^ (pi) = p/ 2 ^/®. Since it is not concave, we do not need to adopt a 
technique to make ^g”^^ (pi) convex as we did previously. We summarize our results in Table 

III. 


j 

77(4) 

''3 

^3 

Po 

j = 1 

p{ 6 p — 2 p^ — ?>) 9 {p — po) 

1 1 

0 0 

1 

0 

Jii « 0.634 

J = 2 

p 2 ^p 2 _ - p)^] 9 (p - Po) 

1 1 

0 0 

1 

0 

1 ^ 0.667 

J = 3 

P® 

p 


2 

21/6 



Table IILSummary of 

and for 

Pi 


B. Case p 2 

In this subsection we would like to quantify the entanglement of p 2 - Above all, we should 
say that Table I implies 

d*’ (P2) = eh (P2) = d"" (P2) = eP (p 2) = o, (2.15) 

because p 2 = p|<h 2 )(*h 2 | + (l —p)|W 4 )(W 4 | itself is an optimal decomposition for those entan¬ 
glement measures. This fact is due to the fact that and cannot detect both |<h 2 ) 
and IW 4 ). 

Let us now compute {p 2 ) and {p 2 )- It is straightforward to show 

[Z 2 (p, if)] = - 2 \/ 6 (l - (2.16) 
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FIG. 2: (Color online) Plot of the p dependence of (a) [Z 2 {p.,p)] in Eq. (I2.16P and (b) 

[Z^{p, ip)] in Eq. (I2.29p . We have chosen (p from 0 to 27r as an interval 0.1. The thick dashed 
lines correspond to F^f^ {p 2 ) in Eq. (I2.25p and (ps) in Eq. (I2.34D . These figures indicate that 
Eq. (12.251) and Eq. (12.3411 are convex hull of the minimum of the characteristic curves. 

where \Z 2 {p.,(f)) is given in Eq. (12.ip . We notice that F^f'^ [Z 2 {p,ip)] has a nontrivial zero 

{p = 0 ) 

_ (2\/6) / _ 

l + (2x/6)2/3 ^ ^ 

Thus, F^'^ (P 2 ) vanishes at 0 < p < po because one can hne the optimal decomposit 


ion 


ft(p) = -P2(P0) + ( 1 - - ) |W4){W4|, 

To V To/ 


(2.18) 


where 


P2(to) = g 


IZ2 (po,0))(^2 (To, 0 ) 1 + 1^2 (po, Y^)(^2 (to, 1 + 1^2 (^Po, Y^)(^2 (^To, y 

(2.19)' 


As the previous cases, we adopt, as a trial, the optimal decomposition at po < p < 1 as 


P2(t) = g 


| Z 2 ( p , 0 ))( Z 2 ( p , 0 ) 1 + 1^2 ( p , y 1)(^2 ( t ) y ] 1 + 1^2 ^ p , y ))(^2 Y 


( 2 . 20 ) 


Then F^'^ (P 2 ) becomes gi{p), where 


9i{p) = -p^^^ P^'^^ — 2\/6(l — p)^'^^ 
9 L 


( 2 . 21 ) 


However, p/(p) is not convex at the region p > p* ~ 0.9196. Thus, we should adopt the 
technique previously used again to make p/(p) convex at the large-p region. 
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Now, we define pi such as po < pi < p*. The parameter pi will be determined later. At 
the region pi < p < 1 we adopt the optimal decomposition for [p 2 ) as a following form: 

P2(p) = T^I*2><*2l+-r^'>2(Pi). (2.22) 

1 - Pi 1 - Pi 

where 


P2(Pl) = 2 


2ti 


27r 


Za (pi,0))(Z2 (pi,0) l + lZa (pi,— ))(Z2 ( Pi, — ) 1 + 1^2 ( Pi, — ))(Z2 ( Pb yj , 

(2.23) 


dvr 


47r\ 


Eq. fl2.22p leads F'f'’ (pa) to gn{p) at the large-p region, where 


9ii{p) = 


P-Pi ^ 1 - p 


1 - Pi 1 - Pi 


p^ - 2\/6pf^(l -pi)^/^ 


(2.24) 


As expected gii{p) is convex at pi < p < 1. The parameter pi is determined by = 0, 
which yields pi ~ 0.86T^. Thus, (pa) can be summarized as 


h-' (P2) 


0 

' gAp) 
gii{p) 


0 < p < Po 
Po < P < Pi 

Pi < p < 1. 


(2.25) 


In order to conhrm again that Eq. fl2.25|) is correct, we plot the p-dependence of the 
characteristic curves (solid lines) in Fig. 2(a) for various p. Our result (12.251) is plotted as a 
thick dashed line. This hgure shows that our result fl2.25p is a convex characteristic curve, 
which strongly supports that our result fl2.25p is correct. 

Now, let us compute (pa). At 0 < p < po, (pa) should be zero due to Eq. (12.181) . 
If we adopt Eq. (12.201) as an optimal decomposition (pa) = gY^ip) is obtained. However, 
it is not convex in the full range. Therefore, we have to choose 


P2(p) = Y ^I^2)(<h2| + -^- -p2{Po) 

1 - Po 1 - Po 

as an optimal decomposition, which results in 

(P2) = 0{p-po) Q 


1/3 


P-Po 
1 -Po' 


(2.26) 


(2.27) 


The parameter pi is obtained by an equation 6pi(4pi — 3)^ = (1 — Pi)(l + 2pi)^. 
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C. Case p 3 


In this subsection we will compute the entanglement of ps = p|<h 3 )($ 3 | + (l — p)|W 4 )(W 4 |. 
Since and cannot detect both |<h 3 ) and IW 4 ), it is easy to show 


d"’ (ft) = sf (ft) = d*’ (ft) = ed (ft) = 0- 


(2.28) 


Now, let us compute (ps) and (ps). For \Zs{p,(p)) in Eq. fl2.ip it is possible to 
show that [Z^{p, p)] reduces to 


[^3(P, p)] =P^ p)e^"^ 


Eq. fl2.29p implies that [■Z^ 3 (p, p)] has a nontrivial zero {p = 0) 


(2.29) 


Po = p = 0.6. 
5 


(2.30) 


Thus, (P 3 ) should be zero at the region 0 < p < po and its optimal decomposition is 


P3(p) = -psipo) + |W4)(W4|, (2.31) 

Po \ Po / 

where 

PsiPo) = ^ [|^3(Po,0))(^3(Po,0)| + |^3(Po,7r))(Z3(po,7r)|]. (2.32) 

If we adopt the optimal decomposition at po < p < 1 as a form 

Psip) = ^ [|^3(p,0))(^3(p,0)| + |^3(p,7r))(Z3(p,7r)|] , (2.33) 

the resulting (P 3 ) becomes |p^ (p — |). Since this is convex, we conclude 

-^ 3 ^^ (P 3 ) = 0{p - Po)^p® (^ ~ ^) ■ 

In order to prove that Eq. fl2.34p is correct we plot again the characteristic curves (solid 
lines) and our result fl2.34p (thick dashed line) in Fig. 2(b), which supports that Eq. fl2.34p 
is optimal one. 

Finally, let us compute (ps). If we take Eq. (I2.33p as an optimal decomposition for 
^ 3 ^^ (Ps) at Po < p < 1, the result is not convex in the full range of this region. Thus, we 
should choose 

Psip) = j —^|*^> 3 )(*^> 3 | + 7 — -PsiPo) (2.35) 

1 - Po 1 - Po 

as an optimal decomposition, which simply results in the right-hand side of Eq. fl2.1ip with 
Po = 3/5. 
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III. DISCUSSION AND CONCLUSIONS 


We compute the three-kinds of true four-way entanglement measures {j = 1,2,3) 
and their corresponding linear entanglement monotones {j = 1,2,3) analytically for 
four-qubit rank -2 mixed states pj = p|<hj)(<l)j|-|-(l — p)|W 4 )(W 4 | (j = 1,2,3). All optimal 
decompositions consist of 2, 3, 4, and 5 vectors. 

Our results can be used to hnd many different mixed states, which have vanishing en¬ 
tanglement. For example, let us consider with po in Eq. fl2.17l) . Let us represent, for 
simplicity, |<F 2 ) and IW 4 ) as 



Imagine the two-dimensional space spanned by |d> 2 ) and IW 4 ) represented by a Bloch sphere. 
Then, the states in the Bloch sphere can be expressed as p = 1(1 -|- r ■ cr), where |r| = 1 
and |r| < 1 denote the pure and mixed states, respectively. In this representation the Bloch 
vectors of |‘F 2 ), IW 4 ), and \Z 2 {po,^)) are 


r(<I) 2 ) = (0,0,1) r(W 4 ) = (0,0,-l) (3.2) 

r{Z 2 {po,p>)) = (- 2 y/po(l -Po) cosp, -2a/po(1 - Po) sinp, 2 po - !)■ 

Thus, any states located in the tetrahedron , whose vertices are (0,0,—1), 
(-2^Po(l -po),0,2po - 1), (^po(l - Po), -\/3po(l - Po), 2po - 1), and 

(\/Po(l “ Po)i ^ 3 ^ 0(1 ~ Po)) 2po — 1) in fhe Bloch sphere representation, have vanish¬ 
ing and . 



C (concurrence) 

r (three-tangle) 

Pi 

i (1 - 2v^-p) 0(ai-p) (ai = (\/2 - 1)2) 

0 

P2 

-P)) ^(«2-p) («2 = 1) 

? 


Cab = |(l - 2x/p - p) 0(ai - p) 


P3 

Cac = Cad = Cbc = Cbd 

TACD = TbCD = 0 


- ^1 P Vp(2 p)) 9{a^ p) (^03 - 2 j 

tabc = TaBD =? 


CcD = \\l - ^ 1 + ^/pi2 - p) + \Jl- y/p(2-p)^| 



Table IV:Entanglement for sub-states of pj (j = 1, 2, 3). 
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One can use our results to discuss the monogamy properties [3^ of entanglement. For this 
purpose, however, we should compute the entanglement for the sub-states of pj (j = 1, 2, 3). 
The entanglement of the sub-states is summarized at Table IV. As this table shows, some 
three-tangle, at least for us, cannot be computed analytically. This is because still we do 
not have a closed formula for computing the three-tangles. 

As mentioned above, there are nine SLOCC classes in the four-qubit system. Therefore, 
many rank-2 states can be constructed by choosing different classes. If, for example. Gated 
and T 70 J are chosen, one can construct the rank -2 state 

=Pl®j)(®il + (l -P)I0(^I, (3-3) 

where |^) = (|0000)-|-|1011)-|-|1101)-|-|1110))/2. Probably, our calculation procedure enable 
us to compute the entanglement of 7^2 and tts although we have not checked it explicitly. For 
TTi, however, our procedure does not seem to work because of (^|*Fi) 7 ^ 0. Using Table II 
one can construct many higher rank states. If, for example. Gated, and La^ are chosen, 
one can construct the rank-3 state such as 

(^3 =Pl®i)(^il+?IO(^l + (l-p-?)lh)(hl, (3.4) 

where \'q) = (|0001) + |0110) + |1000))/\/3. However, it seems to be highly difficult to 
compute the entanglement of higher rank states. 

The most remarkable achievement and novelty of this paper is deriving the entanglement 
of four-qubit mixed states using an analytical approach. Thus, our result may serve as a 
quantitative reference for future studies of entanglement in quadripartite and/or multipartite 
mixed states. 
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